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Abstract 

Continuously symmetric solutions of the Adler-Bobenko-Suris class of discrete integrable equa- 
tions are presented. Initially defined by their invariance under the action of both of the extended 
three point generalized symmetries admitted by the corresponding equations, these solutions are 
shown to be determined by an integrable system of partial differential equations. The connection 
^ ' of this system to the Nijhoff-Hone-Joshi "generating partial differential equations" is established 

and an auto-Backlund transformation and a La:x pair for it are constructed. Applied to the HI 
and Ql members of the Adler-Bobenko-Suris family, the method of continuously symmetric 
reductions yields explicit solutions determined by the Painleve trancendents. 

T— I ■ 

^ ■ 1 Introduction 

\^ • The study of integrable discrete systems has a long history going back to work in the late seventies and 

I early eighties [ll[2l[3l[4j. At this point, it is acknowledged that, most of the well known integrable discrete 

^SJ ■ systems are characterized by their "multidimensional consistency". This means that, such a system 

' may be imposed in a consistent way in a multidimensional space. This property seems to incorporate 

automatically two integrability aspects of this kind of systems, in the following sense: Multidimensional 
consistency allows one to derive algorithmically a Backlund transformation, as well as, a Lax pair for 
the difference equations under consideration O [6l [7] • 

Recently, Adler, Bobenko and Suris (ABS) classified the scalar lattice equations which are multidi- 
mensionally consistent and possess the symmetries of the square and the tetrahedron property, as well 
' [8] . Subsequently, they classified the lattice equations having the consistency property in a more general 

framework [3]. 

The equations covered by the ABS classification [8' have already attracted the interest of many in- 
vestigators and several results pertaining to them have already been published, including exact solutions 
[lOl [TT] , Backlund transformations [12] , symmetries [131 [HI US [El a-^d conservation laws [T^ . 

In this paper, we focus on the symmetry properties of the ABS equations and show how a particular 
class of reductions provide a natural interplay between them and certain non-autonomous systems of 
partial differential equations. The means to explore this link is provided by the pair of extended three- 
point generalized symmetries admitted by the equations of the ABS class [TS". 

More specifically, we study the continuously invariant solutions of the systems under consideration. 
We use the term "continuously invariant solutions" for the solutions that remain invariant under the 
action of both of the extended three-point generalized symmetries admitted by the corresponding equa- 
tion. We show that these solutions are determined by a system of differential-difference equations, 
which involves six values of the unknown function, u. The elimination of three of these values leads to 
an equivalent system of partial differential equations, which involves the remaining values of the 

dependent variable. 
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Among the other advantages offered by the general framework; of the continuously invariant solutions 
developed in this paper is the fact that it allows us to derive easily some of the integrability properties 
of In particular, it enables us to construct an auto-Backlund transformation for this system, as 

well as, a Lax pair. 

The implementation of this general framework to the equations HI and Ql of the ABS family 
leads to explicit solutions, constructed using symmetry reductions of the corresponding S[u] systems. 
These solutions are determined by quadratures from the continuous Painleve V and VI equations, but 
may also be regarded as resulting from reductions which lead to discrete Painleve equations [131 118j . 
In this fashion, a new connection between discrete and continuous versions of the Painleve equations is 
revealed. 

Another important aspect of system is that it leads to what has been termed as generating partial 
differential equations. The archetypical example of such equations is the regular partial differential 
equation (RPDE), introduced by Nijhoff, Hone and Joshi in [19]. These authors showed that the RPDE, 
which encodes the entire hierarchy of the Korteweg - de Vries (KdV) equation, is related to equation 
HI of the ABS family. In the present paper the above result is rederived, but by a completely different 
method, which also allows its immediate generalization. Specifically, wc show that, not only HI, but 
also H2, H3 and Ql are related to the RPDE, and establish this relation in a systematic fashion, using 
the properties of the corresponding 

The present paper is organized as follows. In Section 2, we first introduce the notation used in the 
sections that follow. Then, we present the main characteristics of a wider class of lattice equations 
containing all the members of the ABS family, along with an auto-Backlund transformation, B^j, for 
each member of the latter. 

Section 3 deals with the solutions of the equations of ABS family which remain invariant under the 
action of the two extended three point generalized symmetries admitted by the above equations. These 
solutions are determined by a system of differential-difference equations which we prove to be equivalent 
to the integrable system In the same section we prove that, the class of continuously invariant 

solutions is closed under the Backlund transformation B^. Exploiting this result, we derive two items 
revealing the integrability of system itself, namely an auto-Backlund transformation and a Lax 
pair. 

Sections 4 and 5 contain the application of the general results of Section 3 to the ABS equations HI 
and Ql s^o- Specifically, we construct symmetry reductions of the corresponding T,[u] systems, in terms 
of which, explicit solutions of the above equations are determined. 

Section 6 deals with generating partial differential equations and the detailed analysis of system 
corresponding to equations H1-H3 and Ql is presented. In particular, we show that systems T,[u] for 
HI, H2 and Ql are related, through a contact transformation, to RPDE. Also, we derive the connection 
of S[u] for H3 to RPDE. Finally, we relate our results to the ones of pl^, where the connection of HI, 
H3s=o and Qls=o to RPDE was presented from a different point of view. 

The concluding section contains an overall evaluation of the presented results and various perspec- 
tives. 

2 Notation and the Adler-Bobenko-Suris equations 

We first introduce the notation that will be used in what follows. In addition, we present those properties 
of the ABS equations that will be used in the next sections. 

A partial difference equation is a functional relation among the values of a function u : Z x Z — > C 
at various points of the lattice, which may also involve the independent variables n, m and the lattice 
spacings a, (3, see Figure [TJ i.e. a relation of the form 

(f(n, m, U(o^o),U(i,o),"(o,i)' ■ ■ ■ J"'/') = 0- (1) 
In this relation, denotes the value of the function u at the lattice point [n + i,m + j), e.g. 

Uio.fi) = u{n,m), U(i_o) = u{n+l,m), U(o,i) = u{n,m+l), 
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and this is the notation that wiU be used for the values of the function u from now on. 



"(0,1 



"(1,1) 




"(0,0) "(1,0) 

Figure 1: An elementary quadrilateral 

The analysis of such equations is facilitated by the introduction of two translation operators acting 
on functions on Z^, defined by 

(.yS'^'^u) = U(fe,o) , CyS^'^u) = U(o,fc) , where keZ, 

V / (0,0) V / (0,0) 

respectively. 

The equations of the ABS family belong to a wider class which contains all the equations of the form 

(5("(o,o),"(i,o),"(o,i),"(i,i);a,/3) = 0, (2) 
where the function Q satisfies the following requirements: 

1. It does not depend explicitly on the discrete variables n, m. 

2. It is affine linear and depends explicitly on the four values of the unknown function u, i.e. 

l9u(.,,)Q(u(0,0)>"(l,0):"(0,l):"(l,l);a,/5) 7^ 

and 

^«(i.,)Q("(o,o), "(1,0), "(0,1): "(1,1) = 0, 

where i, j — 0, 1. 

3. It possesses the symmetries of the square (D4-symmetry): 

'3("(o,o),"(i,o),"(o,i),"(i,i);a,/3) = e(3(u(o,o):"(o,i):"(i,o)>"(i,i);/5,a) 

= crQ("(i,o),"(o,o),"(i,i),"(o,i);a,/3) , 

where e — ±1 and a = ±1. 

The affine linearity of Q implies that one can define six different polynomials in terms of the func- 
tion Q [SI ISl US] J four of them assigned to the edges and the rest to the diagonals of the elementary 
quadrilateral where the equation is defined, see Figure [5) 

A polynomial hij assigned to an edge or a diagonal depends on the values of u assigned to the 
end-points of the corresponding edge or diagonal, respectively, as illustrated in Figure [21 and is defined 

by 

hij — i ■ Q Q ,ij Q ,i Qj 7 ^ ~/~ 3 1 ^} j f ; ■ • ■ 1 4, 

where Qa denotes the derivative of Q with respect to its i-th argument and Q ij the second order 
derivative Q with respect to its i-th and j-th argument. The polynomials hij are quadratic in each one 
of their arguments. Moreover, the relations 



hi2 h^i — his — hi4 h 



23 



(3) 



hold in view of the condition Q = 0. 

On the other hand, the symmetries of the square imply the following: 
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/l24 
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W(0,0) U(1,0) 



Figure 2: The elementary quadrilateral and the polynomials 



1) The polynomials on the edges have to be of the form 



hi] 



h{x,y;a,P), \i-j\ = 1 
h{x,y;P,a), \i-j\ = 2 



(4) 



where h is quadratic and symmetric in its first two arguments. 

2) The two diagonal polynomials are identical and 

hiA = h23 = G{x,y;a,f3), (5) 
where G is quadratic, symmetric in its first two arguments and symmetric in the parameters. 

2.1 The Adler— Bobenko— Suris equations 

In order to make our presentation self-contained, we first list all the members of the ABS classification 
and the notation that we will use in the next sections: 



HI (w(o,o) - "(i,!)) (^(1,0) - W(o,i)) - a + (3 = 

H2 (u(o,o) - U(i,i))("(i,o) - "(0,1)) + (/? - Q;)('«(o,o) + "(i.c) + "(o.i) + '"(i,!)) 
-a^ + /32 = 



(6) 



(7) 



H3 a{u(Qfi)U(ifi) + U(o,i)M(i,i)) - /?(w(o,o)W(o,i) + W(i,o)M(i,i)) + 5 (a^ - 0^) = 



(8) 



Ql a(w(o,o) - ■"(o,i))(w(i,o) - "(i,!)) - P{u(o,Q) - ■"(i,o))(w(o,i) - u^i,i)) + S'^a/3{a - /3) = (9) 



Q2 a(u(o,o) - U(o,i))(m(i,o) - ^ /5("(o,o) - W(i,o))(u(o,i) - "(i.i)) 

+a(3{a - /?)(m(o,o) + U(i,o) + M(o,i) + W(i,i)) - "/^(Q! - - a/3 + /3^) = 



(10) 



(35 {P^ - Q!^)(«(o,o)M(i,i) + W(i,o)W(o,i)) + - l)(w(o,o)W(i,o) + U(o,i)"(i,i)) 



-a(/3^ - l)(u(o,o)W(o,i) +M(i,o)W(i,i)) 



J^(a^-/3^)(a^-l)(/3^-l) 
4q/3 



<5^ aoU(o,o)W(i,o)W(o,i)'"(i,i) 

+ai(u(o,o)U(i,o)W(o,i) + "(i,o)"(o,i)"(i,i) + "(o,i)'"(i,i)'"(o.n) + W(i,i)M(o,o)W(i,o)) 
+a2(u(o,o)W(i,i) + W(i,o)W(o,i)) + a2(M(o,o)W(i,o) + W(o,i)"(i,i)) 
+a2(M(o,o)U(o,i) + W(i,o)W(i,i)) + 03(^(0,0) + W(i,o) + W(o,i) + W(i,i)) +ai = 



(11) 



(12) 
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The Oi's appearing in the last equation are determined by the relations 

flo = a + 6 , ai = —a/3 — ba , a^— a0^ + ba^ . 
ab{a + b) 2 (n 2 52 \, _ ab{a + b) 



a2 



ap -[2a a2 = — -- + 6a -(2/3 - -J a , 



2(a-/3) V 47 ' 2{(3 - a) 

93 92 gl 

«3 = Trao - -rOi , 04 = — oo - 3301 , 
2 4 lb 

where 

= p(q!) , 6^ = p(/3) , = 4x^ - 322; - ffs ■ 

The main characteristic of all of the above equations is their integrability, which is understood as 
their being multidimensionally consistent. From this property it follows that [8]: 

i) The polynomial related to the edges, /i, can be written as 

h{x,y;a,f3) = k{a, f3) f{x,y,a) , 
where the function k(a, (3) is antisymmetric, i.e. 

k{(3,a) = -fc(a,/3). 



ii) The discriminant 

is independent of the parameters a, (3. 

Ill) The functions /, G and k determining the polynomials hij can be specified explicitly and, for 
convenience, are given inlAl 

To the above properties of the ABS equations, one can add the following two [7], which will also be 
used in the symmetry analysis to be presented in the following sections. 

iv) They define their own auto-Backlund transformation. The latter is specified by the following rela- 
tions 

B,(«,u,A) ( Q("(o.o),"(i.o),^(o.o),ii(i,o);a,A) = _ 
i 'y(w(o,o),"(o,i)>'"(o,o),'«(o,i);P'''^) = 

v) If {vP ^v} ^v? ^u^'^} is a quartet of solutions related by the Backlund transformation B^, then their 
superposition (Bianchi diagram) is expressed by the condition 



3 Symmetry reductions 

In this section we present the general framework of particular symmetry reductions of the ABS equations. 
More specifically, we study solutions of these equations which remain invariant under the action of both 
of the extended three point generalized symmetry generators, under the assumption that the unknown 
function depends continuously on the lattice parameters a, /3. 

We first show that, invariant solutions of the above kind are determined by a system of differential- 
difference equations. The latter turns out to be equivalent to an integrable system of partial differential 
equations, S[u]. The integrability of S[u] is established by the construction of its auto-Backlund trans- 
formation, Be. This transformation provides the means for deriving a Lax pair for system as well. 
These integrability aspects are the subject of the second part of this section. 
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3.1 Continuous symmetry reductions and system : general considera- 
tions 

Let us recall that [13] , every integrable lattice equation 

^(^(co), "(1,0), '«(o,i),U(i, !);«,/?) = (14) 
admits a pair of three point generalized symmetries generated, respectively, by the vector fields 



^("(0,0) , U{i,o) , "(-1,0) , a) ^«(o,o) ' 

^("(0,0)7 "(0,1), "(0,-1), 9u^o,o) ■ 



(15i) 
(15ii) 



It also admits a pair of extended generalized symmetries with respective generators the vector fields 



where 



vi = A(n)i?(M(o,o),"(i,o),"(-i,o),a)9„(„ ,j) + {A{n) - A{n + 1)) r{a) da , 
V2 = B(TO)i?(M(o,o), "(0,1), "(0,-1), /3)9„,Q_o) + (B{m) - B{m + l))r{(3)dp , 

R(u,x,y,K) := -f^^(u,x,K) = + -f^y(u,y,K), 

X — y 2 X — y 2 



(16i) 
(16ii) 

(17) 



A{n), B{'m) are arbitrary non-constant functions of their arguments, and r{x) depends on the particular 
equation under consideration, as specified in the following table: 



Equation 


HI 


H2 


H3 


Ql 


Q2 


Q3 




r{x) 


1 


1 


X 

2 


1 


1 


X 

2 


-L[Ax^-g^x-g^)y^ 



The solutions of (fT4|) that remain invariant under the action of both of the symmetry generators Vi 
and V2 must satisfy the invariant surface conditions 



du 



(0,0) 



da 
du 



(0,0) 



dp 



K{n,a) i?(u(o,o),"(i,o),"(-i,o),a) , 
i(m, 13) i?(u(o,o), "(0,1), "(0,-1), /3) , 



where 



K{n, a) 
L{m,P) 



A(n)-A(n+1) r{a) 
B(rn) 1 



S(m)-B(m+1) r{P) ' 

The compatibility of equations and (fT5)) is expressed by the conditions 

Da (Q("(o,o),"(i,o),"(o,i),"(i,i);a,/3)) = 0, 
D/3 (Q("(o,o), "(1,0), "(0,1), "(1,1); = 0, 
dp (9qU(o^o)) = da (5/3U(o,o)) , 

where and denote the total derivative operators with respect to a and respectively, i.e 

1 Q,. 1 

Dq := da 



(18i) 
(18u) 

(19) 



(20i) 
(20u) 
(20in) 



E 



du 



da 



i,J=0 



Written out explicitly, by using the expressions for Q a and Q,/3 following from the determining 
equations for the symmetry generators vi and V2, respectively, conditions (|20i|) . (|20iip imply that A{n) 
and B{m) must be affine linear. Without loss of generality, we choose them to read as follows: 



A(n) = n , B(m) 



m . 
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Condition (|20iii|l . on the other hand, imposes no further restrictions, because it holds identicaUy. 
This follows from the fact that, the commutator of the two symmetry generators v„, produces a 
trivial generalized symmetry [7]. 

Thus, the solutions of the ABS equations which are invariant under the action of both vi and V2 
are determined by the differential-difference system 

<3("(o,o)7"(i,o):'"(o,i)'"(i,i);"'/') = 0' 
— + ni?(w(o,o),M(i,o),U(-i,o),a) = 0, (21ii) 

r{P) + mi?(M(o,o);M(o,i)' "(0-1)^/3) = 0- (21iii) 

These solutions will be referred to as continuously invariant solutions. 

System (j21[) involves the values of the unknown function u at six different points of the lattice. 
One could eliminate any three of these values and get an equivalent system of partial differential 
equations involving the remaining ones. We choose to eliminate the values U(_i^o)i W(o,-i) and W(i,i), 
and this leads to the following result. 

Proposition 3.1 

Every continuous invariant solution is determined by the system of partial differential equations 
9u(i,o) _ G{u{i,o),U(o,i)) 9u(o,o) , "^/(u(o,o),W(04)'/3)^ f G(m(i,o), U(o,i)) 

-Or, 



d(3 A:(a,/3)/(u(o,o),U(o,i)'/^) 2r(/3)fc(a, /?) V/(mo).^i(o,i),/3); ' 

"^«(l,0) 



da fc(/9,Q;)/(M(o,o),U(i,o),a) da 2r{a)k{/3,a) V/(«(o,o), M(i,o) , a 



5^"(o,o) _ . 9-U(o,o) 9u(o,o) , /(w(o,o),"(i,o),a) ( 9-U(o,o) , n 9u(o^o) 



dadji da 8/3 ' 2fc(a,/3) \r{p)"^ da ' r{a)' ^' dp 

nm/(u(o,o),U(i,o):a) 



4r(a)r(/3)fc(a,/3) 



where 



■ _ i' /,f(o.o)('»(o,o)^"(i,o)^Q^) /('»(o,o):'"(i,o)>Q^) a f G(m(i,o),U(o,i 



/(w(o,o),W(i,o),a) fc(a,/?)/(u(o,o),W(o,i)'/3) \/(w(o,o), "(1,0): ") 

/,«(0,i) (^(0,0) J ^f(0,l) I /3) G(u(i,o) . ^^(0,1)) ^,"(0,1) ('"(1,0) > ■"(0,1) ) 



0) 



/("(o,o),U(o,i),/3) /('i(o,o),U(i,o),a) /("(o,o), U(i,o), ") 



A i( a\p, l^ ft \\ /("(0,0),M(l,0),aK /<^,«(i,o)("(l,0):"(0,l))\ 

^3 = fc a, a„ (In / W(o,o) , M(i,o) , a ) - Tr^^-^-^i^^d o) — FT^ ^ 

' ' / W(o,o),M(o,i),/3 ' V / W(o,o),'«(i,o),a / 



G(M(1,o),U(0,1)) „2 /, ^^ 

7?— — — — 5Y^«(i,o) (ln/(M(o,o),U(i,o),a)) , 

/(,W(0,0): ^'(0,1): PJ 
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and 



. _ c, „ //,«(i.o)("(o,o),'«(i,o),a)G(u(i_o),W(o,i)) 

^4 - I f7;;~~";;~~;vi lz,„ (,u(i,o),w(o,i)J 

V /l"(0,0)j^^(l,0)7«j 

, « /, x;, <^,"(1.0)("(1,0)>'»(0,1)) 

+ d„ (ln/(u(o,o),M(o,i)'P)j /("(o.o),'«(i.o),ajc'«(io) — FT ^ 

\ / 1^^(0,0) J '"(1,0) I «j 

- ^«(o,i) (ln/(M(o,o),W(o,i),/3)) 5^(10) (ln/(u(o,o),U(i,o),a)) G(u(i,o) , U(o,i)) • 

The above system, which will be denoted by S (u(o,o) j ""(1.0)1 '^'(o.i) ! Qf; '^i or, simply is symmet- 

ric: 

S (u(o,o),'"(o,i),U(i,o);^,a;"i,n) = E (u(o,o), W(i,o) , "(0,1); a, ^; ■ 

Proof 

The first equation of system results by eliminating the values U(o,_i) and M(i__i) from equation 
((2T111I) and 

Q(M(o,-i):"(i,-i),'"(o,o),W(i,o);a,/3) = 0- 
Using the affine linearity and the symmetries of Q, the last equation can be written as 

Q,«(0.i)M(0,-l) +Q 
"(1,-1) : l^^J 

V, U(o, 1)41(1. i)'"(0,-l) + V,U(i.i) 

where the arguments of Q(w(o,o): "(1,0): "(0.1)1 ""(Li); Q!? /3) have been omitted and Q and its derivatives 
are understood to be evaluated at ""(0,1) — "(1,1) = 0- 

We now solve (j21iiip and its shift in the n direction with respect to W(o^_i) and W(i respectively, 
and substitute the results into equation (|22p . The resulting equation, combined with the relations 

^2 _ /('»(o,o),"(i,o),")G'(m(i,o),M(o,i)) 

'"<^'^' " /("(i,o),"(i,i),/3) 
„ „2 _ fc(Q,/3)/(u(o,o), "(1,0)7 «) (^,^0 1) ("(1.0), "(0,1)) - /,«(i 1) (m(i.o),M(i,i),/3)) 

" /("(1,0),"(1,1),/?) 

which hold in view of Q 0, yields the first member of 

The second equation of T,[u\ results in a similar manner. It is also easily verified that, the first two 
equations of are symmetric, i.e. the one is mapped to the other under interchanges 

"(1,0) < — '"(0,1), a< — > /3, n< — > m . (23) 

In order to obtain the third member of one first solves the second equation of and its 

shift in the n direction for 9/3U(i_o) and c}^U(_i^o)- One then substitutes the result into the derivative of 
equation (|21iip with respect to /3. From the resulting equation, one arrives at the third member of T,[u] 
by using the expressions for U(i.o) — "(-i.o)> G(u(_i 0) , U(-o_i)) and its derivatives provided by equation 
(|21iip and the relatiorQ 

^ . ^ . ^ ^ g("(l,0),"(0,l)) -g("(-l,0),"(0,l)) 

^,"(1 o)("(l,0),"(0,l)j + (^,U(_i o)("(-l,0),"(0,l)j - ^ , 

"(1,0) — "(-1,0) 

and its differential consequences, respectively. 

Finally, differentiating equation l|21iiip with respect to a and following an analogous procedure, one 
arrives at an expression which is identical to the third member of under the mapping (|23p . □ 



^This relation holds identically, i.e. without taking into account the equation Q = 0, thus we can differentiate it 
assuming that the corresponding values of u are independent. 
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As already noted, one may choose to eliminate any other triad of the values of u involved in equations 
(|2ip . In this fashion, compatible systems of partial differential equations can be constructed involving 
the triplets (^(o^o), "(-i.o), "(o.i))' ("(o,o), "(i.o), W(o,-i)) and (u(o^o), "(-1,0), "(o-i)), respectively. It 
turns out that, the corresponding systems are given by 

S(u(o,o),W(-i,o):M(o,i);a:/3; -n,m) , I](u(o_o): W(i,o), M(o,-i); a, /3; -m) 

and 

5^(^^(0,0): "(-1,0), W(o,-i); -n, -m) , 

respectively. 

System I](u(o,o)7 ''^(i.o), "(0,1)! Q^, P] '^j ™) and the last three are compatible, in the following sense, cf. 
Figure[3l If we eliminate 9^M(i_o) (respectively 9/3it(_i 0)) from systems I](ri, m) and — m) (respec- 
tively I](— n, m) and !](—«, — m)), then we will end up with (I21iii|) . On the other hand, the elimination of 
9aU(o,i) and 9qU(o.-i) from systems I](n, m), S(— n,m) and S(n, — m), I](— n, — m), respectively, leads 
to (|21ii|) . Finally, the elimination of 9a9/3U(o,o) from any two of the four E[it]'s mentioned above results 
in (I2TTI1 . 

"(0,1) 

E(— n,m) E(n, m) 
U(_i,o) • •"(1,0) 

"(0,0) 

E(— n, — m) E(ri, — to) 

"(0,-1) 

Figure 3: The values of u and the compatible systems S 



3.2 Integrability of system I][u] 

We have already characterized system Efu] as integrable. To support this characterization, in the present 
subsection, we construct an auto-Backlund transformation and a Lax pair for the above system. 
To this end, let it first be noted that the fact that every integrable lattice equation 

'9("(o,o),"(i,o),"(o,i),"(i,i);a,/3) = (24) 

admits generalized symmetries and extended generalized symmetries with generators the vector fields 
given in (jl5i[) and (jl6ip respectively, has the following consequences. 

1. The vector fields 

v„ = -R ("(0,0), "(1,0), "(-1,0), a) 9u^o,o) + ("(0,0), "(1,0), "(-1,0), 5«(o,o) 

and 

vi = A{n) v„ + {A{n) - A{n + 1)) r(a) , 
are symmetry generators of the first of the equations making up the auto-Backlund transformation 
Bd, i.e. of (3(m(o,o), "(1,0), "(0,0), "(1,0)! Q!, ^) — 0- Therefore, relations 

Vn^ (Q("(o,o),"(i,o),"(o,o),"(i,o);a, A)) = 0, 
v^^^ (Q("(o,o),"(i,o),"(o,o),"(i,o);a, A)) = 

hold in view of Q(u(o,o), "(1,0), "(0,0), "(1,0); A) = 0. 
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2. The vector fields 

Vm = ("(0,0) J ""(0,1) '""(0,-1) 7/5) C^«(o.o) + (^'(0,0)7'"(0,1)'"(0,-1)'^) ^"(0,0) ' 

V2 = S(m) v,„ + (B(m) - B(m + 1)) r(/3) , 

are symmetry generators of the second of the equations of the auto-Backlund transformation. As a 
result, the pair of relations 

v^^ (Q(w(o,o),W(o,i)'"(o,o),U(o,i);/3'^)) = 0' 
v^^^ (Q(m(o,o),W(o,i),U(o,o),U(o,i);/^, A)) = 

hold in view of Q{ui^o,o),U{o,i),U{o,o),U{osy-, PA) = 0. 
Using the above observations, one may prove the following proposition. 
Proposition 3.2 

The auto-Bdcklund transformation Bd(u,'U, A) maps a continuously invariant solution u to another 
solution u of the same kind. 

Proof 

It is given in the |B] □ 

An immediate consequence of the this result is described in the following proposition. 
Proposition 3.3 

If u is a continuously invariant solution, then the system 

Q(w(o.o),W(i,o):'"(o,o)7'«(i,o);", = 0, (25i) 
Q(u(o,o),U(o,i),U(o,o),'«(o,i);/3, A) = 0, (25ii) 



9u(o.o) _ 1 / i9it(o,o) , f / G(M(i,o),'"(o,o),a, A) 



/ 19 "(0,0) ^ n ^ \ 
+ ^/.«a,o,(mo),"(i,o),«)j 



da k{a,\) \ da 2r{a) ' ^ ' J /(u(o,o), U(i,o), 

- 2k{a,X)r{a) ^'"<^''" ("(i.o), ^(o,o), «, A) , (25iii) 

5u(o,o) _ 1 / 9w(o,o) , _f / G'(u(o,i),W(o,o),/?, A) 



+ ^7(;^/.«(o.i)("(0,0),U(0,l),^) 



8(3 fc(/3,A)V 5/? ■ 2r(/3)^'"<°'^'^-^^^'"^'™'-7 /Ko,o),"(oa),/3) 

in 

- 2fc(AA)r(/3) ^■"<-' ("(04). ^(0,0), /3, A) , (25iv) 

defines a new solution u of the same kind, and conversely. 

The above system, which will be denoted as Bc(u,?i, A), is symmetric 

Bc(u,M, A) = Bc(m, u, A) , 

and defines an auto-Bdcklund transformation of system Y,[u\. 

Proof 

It is given inO □ 
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Remark 3.1 

i) The second pair of equations of system ]Bc(u,{t, A) follows from the two first equations of S[u], via 
the substitutions 

U(i,o) — » ""(0,0) , oi — > A 

and 

'^(O.l) * '"(0,0); P ^ -^J 

respectively. 

ii) The superposition principle of B(j(m, u, A) implies the corresponding one for Bc(u, u, A): 



"(0,0); 


"(0,0);^1''^2) 


= 0, 


"(1,0); 


"(1,0); ^1' "^2^ 


- 0, 


-7/2 

"(0,1); 


"(0,1)5 ^1, 


0. 



Finally, let us consider the following pair of equations 



where 



/("(o,o),"(i,o);a) 



n 



2r{a) 



B 



□ 



$ = yy^ f ^An(„,0)«(0,0) ) <j, ^ (26i) 



. / 5u(o.o) , n \ G(u(i o),'S(o,o),a, A) 

^ (^-^^ + ^7(^/,-(i.0)("(0,0),"(l,0),«) ' 



^,"(1,0) ("(1,0) ; "(0,0) , a, A) , (27i) 



5u(o,o) , ™ f ^ G'(u(o,i),W(o,o),/3, A) 

— aTs o /ON •'."(0,1) (."(0,0); "(0,1); Pj —77 

dp 2r{f3) ) /(m(o,o),"(o,i);/3) 

- ^^:^G„(„^^)(w(o,i),U(o,o);/?;A) , (27ii) 

and A, B and their derivatives are evaluated at W(o,o) = 0- Equations (|26i[) constitute a Lax pair 
for system I][it]. One arrives at this result, essentially, by the inverse of the procedure presented by 
Crampin in [20]. In any case, it can be easily verified directly, by considering each of the ABS equations, 
separately. 

4 Continuous invariant solutions of the discrete potential KdV 
equation 

In the last section, we established the general framework for the special reductions of the ABS equations 
leading to what we called continuously invariant solutions. In the present section, the above results are 
applied to equation HI. The latter will be also referred to as discrete ■potential KdV, in compliance with 
the terminology adopted in pP , cf. also [2l [3] . 
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System E[it] corresponding to HI is made up of the equations 



dui 


Ui - U2 




a- 13 


du2 


Ui - U2 


da 


a- P 




a-P ( 


dadp 





m~ {ui- U2) , (28i) 
n + {ui~U2)-^ \ , (28u) 

where 

U = U(^o.O) , Ul = '"(1,0) , U2 = M(0,1) ■ 

Obviously, the nonhnear system ([28]) is very hard to solve. However, whole families of solutions can 
be obtained, in a systematic way, via symmetry analysis. In what follows, we construct multiparameter 
families of solutions of the above system, using its Lie point symmetries. 

System I^E^ admits a five dimensional group of point symmetries generated by the vector fields [21] 

^1 = da + dp , W2 = ada + Pdfj + udu , 

W3 = a„ , W4 = + du2 , W5 = udu - uidm - U2du2 ■ 

It will be shown that, similarity solutions corresponding to above group of symmetries are determined by 
solutions of the Painleve V and VI equations pi]. For easy reference, we note that, the latter equations 
are given by 



2G + G—l)^ + 



and 



1/1 1 1 



2 \G G-1 G-yJ \y y-l G-y^ 



(30) 



y^y-ir V G2 ■ (G-l)2 {G-yf 



which will be denoted by ^v{y, G{y); a, b, c, O) and ^vi{y, G{y);a, b, c, 5), respectively. 

For the same reason, we list here the symmetry generators of the discrete potential KdV equation 
[T51 115j . which will be used extensively in the rest of this section: 

• Point symmetries 

XI - X2 = (-i)"+"a„,„,„,, X3 ^ (-i)"+™«(o,o)a„,„,„, , 

X4 = ■"(o,o)5«(o,o) + 2a9a + 2f3df3, :x.5 ^ da + df3 . 

• Generalized symmetries 

^(1,0) ~ ^(-1,0) — p) 

" . y, — :^"<«.°) ' ^4 = m V3 - , 

^(0,1) ~ "(0,-1) ^1°; — pj 

V5 = nvi - a^,, vg = mva - 9^. 
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4.1 Solutions related to Painleve V 

We first consider solutions of system (|28p that are invariant under the action of the symmetry generator 
y^i + 2 ^yvz = da + dfi + 2 ^ {udu - uidu^ - U2du2) , ^ G M - {0} , (31) 
where /i may depend on n, m. Such solutions have to satisfy the differential equations 

This implies that, u, ui and U2 are given by 

u(a,P) = Tn^niiy) exp(^z), 

Mi(a,/3) = Tn+i,m{v) exp(-/iz), (32) 
U2{a,[i) ^Tn,Tn+i{y) exp(-/xz), 

where 

y — a — (3 , z = a + P , 

and Tij{y) are arbitrary functions. 

Substitution of these forms into ((28|) yields the following system of ordinary differential equations 

y {K+i.rniy) + ^^Tr^+iMy)} - - {m + j^-^) 

2/ - A*7^«,m+i(2/)} = (n + ^+^)^, (33) 

y^lra^y) - {n + ra)T'^^^{y) = 2^+^-^ - (n-m-My)^r„,™(y), 

where 

— Tn^niiy) ± f^Tn^miy), SS := r„+i,„,(y) - r„^„+i(y), 

and the prime denotes differentiation with respect to y. The analysis of the latter system can be 
summarized as follows. 

Starting from (j33p and using differentiation and elimination, one arrives at a fourth order ordi- 
nary differential equation for T„_m(y), which is omitted because of its length. The order of the latter 
differential equation is reduced by one, using the quadrature 

^ln(T„,„(y)) = M^i^^. (34i) 

dy 1 - Gn^miy) 

Finally, the resulting third order equation can be integrated once to yield 

f2/,G„,„,(y);y,-— ,A,-2AiM ^ (34^^ 

where A is a constant of integration and may depend on the parameters n, m. 

Returning to system ([55]) and using relations , one finds the following expressions for the functions 
Tn+l^m{y), [y)-- 

, . yG'n,m{y) + nGlmiy) - (2« + 1 + K - 2^ly) Gn,m(,y) + n + K + 1 

A^lil - Gn, m[y))Tn,m{y) 

y ^^^^^ ^ - {m + K+ l)Gl^„,{y) + (2m + 1 + At + 2/iy) G„,^(j/) - m ^^^..^ 

"'™ - Gn,rn{y))Gn,m{y)Tn,m{y) 

where k = A/ (2/i). 

In order to satisfy the discrete potential KdV equation, the functions u, ui and U2, given by 
and (|35ip . must be related by shifting appropriately n and m, i.e. 

Ml = yniu) , U2 = ym{u). 
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The above conditions imply that parameter /i must have the form 

M - (-1)"+"t, r e R, 
and functions Tij{y) have to satisfy the foUowing relations: 

r„+l,m(y) = -^n (Tn.miy)) , T'n,m+l(2/) = (T'n,™ (y)) . 

The combination of these conditions with equations ([34]) and (|35i|) leads to the following restriction 
on parameter A: It must be of the form 

\ = P - (-1)"+"t, pen. 

Recapitulating, we can state that, the discrete potential KdV equation admits continuously invariant 
solutions of the form 

uiofi) = r„,„(a - f3) exp [(-1)"+"V x (a + /?)] , r e R , (361) 
where Tn^miv) is given by the quadrature 

|;Mr.,..(„)).(-ir"vi±|^, (3M 

with Gn,m{y) being a solution of the Painleve equation 

(^y, Gn,m{y); Y' - P - (-1)"+"^, -2 j , pen. (36iii) 

Remark 4.1 

The solutions of the discrete potential KdV equation just constructed can also be considered as being 
derived from what is referred to as the asymmetric, alternate discrete Painleve II equation. This ob- 
servation results from the following considerations. By construction, any solution u of the class derived 
above satisfies the differential equation 

do^uf^ofi) + 9/3W(o,o) = 2t(-1)"+™u(o,o), (37) 
as well as, the invariant surface conditions (|21iip . (|21iiip . i.e. 

Ti Tfl 

9aU(ofi) H = 0, 9/3U(o,o) H = 0- (38) 

"(1,0) - "(-1,0) "(0,1) - "(0,-1) 

Elimination of the derivatives of U(o,o) from ([57]) . ([55]) leads to 

77 777 

+ + 2 r (-1)"+" 7.(0,0) = . (39) 

"(1,0) - "(-1,0) "(0,1) - "(0,-1) 

This, however, is the invariant surface condition for solutions of the discrete potential KdV that remain 
invariant under the action of the symmetry generator V2 + V4 + 2tx3 . As shown in [13| , this class of 
group invariant solutions are determined by solutions of the asymmetric, alternate discrete Painleve II. 
□ 

4.2 Solutions related to Painleve VI 

The solutions of system (j28p which remain invariant under the action of the symmetry generator 

W2 + (2/z - l)w5 = adc + (3dp + 2nudu + (1 - 2//) {uidu, + W29„ J , ^ G R - {0} , (40) 
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where n may depend on n and m, must satisfy the differential equations 



a"i,a + f3uij3 = (1 - 2fi)ui , 
aU2,a + /3w2,/3 = (1 - 2^)W2. 



Hence, these invariant solutions must have the form 



(41) 



where 



a 



a l3 . 



leads to the system of ordinary differential 



Substitution of the above expressions into system 
equations 

(1-2/) {2yS'^+^^„,{y) + (2^^ - 1) 5„+i,™(j/) } = 2 {m + ^s/-^} 3S , 
{y~\) {2y5;^„+i(y) - (2/. - 1) 5„,™+i(y) } - 2{nj/ + Vy=e^+^} ^, 
(y - 1) {y'Sl^^{y) + 2/5;^„(2/) - A*'^„,m(2/)} = 2^^^+^^^ + + ny^- , 

where 

■= yS'n^miy) ± f^SnMv)^ ^ — 5'„+i,™(y) - S'„,™+i(?/), 

and the prime denotes differentiation with respect to y. The analysis of the latter system is similar to 
the one described in the previous subsection regarding system (|33l) and leads to the following results. 
The function Sn,m{y) is determined by 



— \n{SnMy)) = FT TT' 

dy y y- Hn,m{y) 



(42i) 



where the function Hn^miy) is a solution of the equation 



'vi y,Hn,m{y)\ 



1 



(42ii) 



In the latter, A stands for a constant of integration which may depend on the parameters n and m. 
The functions Sn+i.m{y), Sn,m+iiy) are given by 



Sn+l,miy) = 



Sn,m+liy) 



y\y - ifKj + 2{2^i - l)y{y - l)i?„.™(i/„,™ - l)i/;.„ 



8fl{2fl- 1) yl/2 (iJ„ ,„ - y) {Hn,ra " 1) S"™,™ 



y^(y - i)^h;,„, + 2(2m - i)z/^(;/ - i)(g»,,n - i)g;,^ + 

8 ^ (2^ - 1) 2/1/2 (iJ„ „j _ y) {Hn^rn " 1) ^f«,m Sn^m 



where we have omitted the argument y of Hn^m and Sn,m- In these relations, summation over the 
repeated index i = 0, ...,4 is understood and the coefficients Ai — Ai{y ,71,771), Bi — Bi{y ,71,771) are 
given by 



AQ{y,7i,7n 
Ai{y,7i,7n 
A2{y,7i,7n 
As,{y,7i,7n 
Ai(y,7i,7n 
Bi{y,n,7n 



2 2 

—TO y , 

y [(m^ - 2A + (n - 2^ + 1)^) y + 2m2] , 

-{71 -2fi+ 1)\2 _ 2 (m^ - 2A + (n - 2/i + 1)^) y + (1 - 2yLf 
2{n - 2/i + l)^y + + 277^ - 2A - (n + 2/i - 1)^ , 
-n(n - 4/i + 2) , 

2/^A4„,,(2/"\TO,n) , i = 0, ...,4. 
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Remark 4.2 

For later purposes, it is noted that, the functions Sn+i,miy)i Sn,m+i{y) may also be considered as being 
determined by the Painleve VI transcendent. Specifically, these functions are determined through the 
quadratures 

_ 1^2m y+H„+i,,„{y) 



3^ In (S'n^m+l(?y)) = 21/ ^ + ' 

where -ff„+i^m(y), Bri,ra+\{y) satisfy the equations 



(44i) 



^vi ^ ^ ,-— ,A,2A^(l-/^)j , (44ii) 

and 

■^Vi U,gwi(y);y,- ^ ^ \ A,2/x(l-A^)j , (44iii) 

respectively. 

This can be proven in the following fashion. Combining relations (j43i|) and (|44ip . we express 
Hn+i.m{y), Hn,m+i{y) in tcrms of Sn,m{y), Gn,miy) and their derivatives. Substituting the result- 
ing expressions into equations (|44iill . (|44iiip . we arrive at (|42iip . □ 

In order to satisfy equation HI, the functions u, ui and U2, determined by (HJ), (jH]) and (|43ip . must 
be such that, ui and U2 result from u by applying the shift operators on u, i.e. 

These conditions imply that parameter ^ must have the form 

A* - i (1 + 2(-l)"+"V) , reM, 
and functions Sij{y) must satisfy the following relations: 

Sn+lMy) = '^niSn.miy)) , S'n,m+l(2/) = (5'„,m(j/)) . 

It is easily verified that, the only consequence of the last conditions is that the parameter A is 
independent of n, m. Hence, we conclude that, HI admits continuously invariant solutions of the form 

a \ , ^^(l + 2(--l)" + "r)/4 



^^(0,0) = ^n,™ ( ^ ) iaPr^^'^'' , r eR, (45i) 



where the function Sn.m{y) is determined by the quadrature 

^ WQ ( 1 + 2(-l)"+'»r y + Hr,,m{y) 

4?/ y-Hn,m[y) 

with Hn^miy) a solution of 

^vi Uh^M y - — 2 ~ 8 2(-l)"+"V)'j , A e R. (45iii) 

Remark 4.3 

Function u, defined in (|45p . satisfies the invariant surface conditions psp and, by construction, the 
differential equation 

(46) 
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Elimination of the derivatives of U(o,o) involved in ([38| and ((46|l leads to 



an 3m 
+ 



"(1,0) ~ ""(-1,0) "(0,1) - "(0,- 

The last equation implies that, solution (|45p is also invariant under the action of the generalized sym- 
metry generator av2 + /3v4 + rx3 . 

Reductions of HI using the above symmetry generator were studied in [13j . while the connection of 
the corresponding similarity solutions to discrete generalized and continuous sixth Painleve equations 
was demonstrated in [18j . On the other hand, similarity solutions corresponding to the symmetry 
generated by av2 + /3v4 were studied in . The latter are contained in the class of solutions given by 
m for T = 0. □ 



5 Continuous invariant solutions of the discrete Schwarzian 
KdV equation 

The Qlg^Q ABS equation, which is also referred to as the discrete Schwarzian KdV [21 , is given by 

a(w(o,o) - w(04))("(i,0) - "(1,1)) - /3("(o,o) - "(i,o))("(o,i) - "(i,i)) = . (47) 

In this section we present continuously invariant solutions of the above equation using similarity solutions 
of the corresponding system 

For this reason, we first write out explicitly and list the algebra of its Lie point symmetries. 
Specifically, S[u] is made up of the equations 



dvi _ VI-V2 "^i"^- "i) (" - "2) - /3 ("1 - "2) 

W ~ a- (3 {v - v^y 

dv2 _ V1-V2 n{v~ vi) {v -V2) + a {vi - V2) || 

da a — j3 (v — vi) 



2 



(48i) 
(48ii) 



d'^v \ ( ( a 3 \ dv dv dv dv \ , , 

+ ""^ ~ "^"5" ' (48ni) 



dadp a — (3\\v~vi v — V2 J da d(3 djd da 



where 



" = "(0,0) , "1 = "(1,0) , "2 = "(0,1) ■ 

Its algebra of Lie point symmetries is three-dimensional and is spanned by the vector fields 



zi = ada + Pdp , Z2 = 9t, + d^^ + dy^ , Z3 = v dy + vi dy^ + V2 d, 



V2 



The solutions of system (j48p remaining invariant under the action of the symmetry generator Zi-f 27Z3 
are determined by the Painleve sixth transcendent. To derive this result, we will make use of the following 
lemma, whose proof is straightforward. 

Lemma 5.1 

The contact transformation 

dv ui — V f , du\ , , 

' ■ ("i-")Tr: , (49i) 



da a \ da 

dv "2 — " / / ^ du\ . 

m + {u2-v)—] , (49u) 



9/3 P \ ' ' d(3, 

vi = ui, W2 = "2 , (49iii) 



maps solutions of system ( fg5|) to solutions of system and conversely. 
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v{a,(3) 




H/2 




= Pn+l,miy) 


-M+1/2 


U2(a,/3) 


= Pn,ra+l{v) Z" 





Thus, we start with the similarity solution of system ([28]) specified by ((4T|) . (|42|) and (|43ip . in which 
we set A = 2^"^, for later convenience. Substitution of this solution to (|49ip and integration of the result 
leads to the following solution of S[w]: 



(50) 



where y = a/ (3, z — a(3, and 

p , N (m + n - 2(£+ ^) + l)(Hn,m{y) - y) Sn+l.mjy) - Hn,m{y)S„,m+l{y) 

'""^^^ " Afl^{Hr,,rniy) - l)^n,™(y) i?„,™(2/) "1 ' ^ ' 

Pn+l,7n{y) = 5„+i,„(y), (51h) 

Pn,rn+l{y) = -Sn^m+i (y) . (51111) 

The functions Pij (y) are determined by solutions of the sixth Painlcvc equation. Regarding Pn+i,m (y) 
and Pn,m+i{y), this property is an obvious consequence of equations (|51ip . 

On the other hand, Pn,m{y) is determined by the solution Hn,m{y) of the Painleve VI equation. 
More specifically, the former is determined by the latter through the quadrature 

UniP.My))-'-^ '-f-^'\ , (52i) 
2y y + Hn,m{y) 

where Hn^m{y) stands for any solution of the Painleve equation 

^vi (y,ff„.„(y); y,-— ,^^-^,2Mi-A*)) ■ (52ii) 



The proof is similar to the one described in Remark HI 
So far, we have shown that, the triad of functions 

v{a,[i) = Pn,ni{y) Z^ , 
Vi{a,l3) Pn+i^r,i{y) Z^ , 
V2{a,/3) = Pn^ra+l{y) Z^ , 

where 

and the Pi^s are determined by solutions of the Painleve VI equation, provide a solution of system 
(|48p . Obviously, all members of this triad are also invariant under the symmetry generator zi + 27Z3, 
i.e. they satisfy the differential equations 

respectively. 

The functions v, vi and V2, as defined above, form a solution of the discrete Schwarzian KdV equation 
provided that 

Vl = yniv) , V2 = ^m(w)- 

These conditions imply that, parameter 7 is independent of n, m, and parameter i must have the form 

£ = ^ {c + n + m + 1) , ceR, 
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in view of which, functions Pij{y) satisfy the following relations: 

Pn+l.miy) = yn{Pn,m{y)) , Pyum+liv) = (Pn.Tniy)) ■ 

As a result, the continuously invariant solutions of Qlg^Q constructed above can be written as 



"^(0,0) 

where 



{a/Sy , 7 G IR, (53i) 



^ln(P„,™(.)) = ^^^tf^, (53ii) 
and Hn^m{y) satisfies the continuous Painlcve VI equation 

^Yi\y,H„^^{y)-—,-—,-{n + m + cf, c G R. (53iii) 

Remark 5.1 

The continuously invariant solution (|53p also satisfies the differential-difference equations 



OOi V(is)) - ^^(-1,0) 

^ ^^(0,0) ^ ^ (^(04) -^(o,o))(t'(o,o) ~t;(o^,i)) ^ ^^^..^ 
W(o,i) - ^(0,-1) 

The latter are but the invariant surface conditions (|21ii|) . (I21iiip . On the other hand, the function W(o.o) 
also satisfies the differential equation 

ai9ai;(o,o) + /3i9/3W(o,o) = 2 7^;(o,o) , (55) 



since it is invariant under the symmetry generator Zi + 27Z3. Using equations (j54ip to replace the 
derivatives of ti(o.o) appearing in the last equation, we conclude that every continuously invariant solution 
must satisfy the following constraint 

(W(1,0) - V(0fi)){v(0fl) - -^(-1,0)) , («(0,1) - V(Ofi)){v{Q.,Q) " ^(0,-1)) „ 

n 1- m = 2 7 7;(o o)- 

w(i,o) - ^(-1,0) w(o,i) - ''^(0-1) 

Reductions of the Schwarzian KdV equation constructed on the basis of this constraint were presented 
in [HI. □ 



6 Continuous invariant solutions and generating equations 

The notion of generating partial differential equations was introduced by Nijhoff, Joshi and Hone in [19j . 
where their archetypical example, the RPDE, was also presented. In the present section, we show that 
system corresponding to several members of the ABS class is intimately related to the above kind 
of equations. Our method of deriving this relation enables us to produce the results of [TH] in a more 
systematic way, as well as to extend these results to other integrable lattice equations. In particular, 
we show that equations H1~H3 and Ql are related to RPDE. 

Remark 6.1 

In order to simplify the resulting expressions, in the present section, we adopt the following notation 
for the corresponding it(i,j) 

U = M(0,0) , Ui = W(i,o) , U2 = W(0,1) ■ 

□ 
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dui 


Ui - U2 


'W 


a — (3 


du2 


Ui - U2 


da 


a — (3 




a-P ( 


dad/S 





6.1 The system S[m] of HI, H2 and Ql 

Let S{U,ui,U2', 6) denote the following system of partial differential equations: 

m-{u^-U2)^^ + ia-P)S'^, (56i) 

n+ iui-U2)^) - {a-/3)S^^, (56ii) 
oa J oa 

One arrives at the above system starting from corresponding to equations HI, H2 and Ql through 
the following contact transformation 

2/(u,iii,Q) U .a -\- n f .iii{u,u-i^a) 
= 2r{a) 

{u,U) := { . (57) 

_ 2 f{u,U2,[i) U_i3 + m {u,U2,0) 
""./S — 2r(/3) 

In particular, ^{u, U) maps 

1. of HI to S'[m;0], 

2. of to ^[t/; S] with (5^ = 1, and 

3. of to S[U;S]. 

In view of these observations, system S[U ; 5] introduced above incorporates the continuously invari- 
ant solutions of the three integrable lattice equations HI, H2 and Ql. What is remarkable is the fact 
that S\U; S\ is also related to RPDE, which has been shown to be a generating equation for the KdV 
hierarchy [19 . 

Indeed, S[U;S\ can be decoupled leading to a fourth order partial differential equation for each of 
the functions involved. To see this, we first solve equation (|56iii[) for the difference ui — U2 to find 

1 // U ad m n \ . 

"'-- = 2(<'-«t757; + i7;-zr:)- '''' 

Substituting the above expression into (j56ip . (|56ii[) . we obtain dfjUi and daU2 in terms of U and its 
derivatives. Then, we differentiate equation (|58p with respect to a and use (|56iip and ((58| to eliminate 
daU2 and ui — U2, respectively. This gives daUi in terms of the derivatives of U. The compatibility 
between the resulting expression and the first equation of S[U;d] leads to the following fourth order 
partial differential equation 

^(a, A -) + 2 - - (2^ - — ^ - - j j = , (59) 



where 



{a,fj,U;n,m) := -L/,aa/3/3 + ^ + TT ^ 71— + 



a -(3 U.p U^a J '"\(3-a C/,„ 
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^'"^ (C/,a + C/,/5 + 2(/5 - a)C/,„/3) 



2(a - /3)3 

f/ 1/11 

2[a-(jyU^l3 2{a-f3) ' \U^a U^p 

When (5 = 0, the last equation reduces to ^(a, /3, U\ n, m) = 0, and this is exactly the equation 
named RPDE [191 HSl- In fact, even when (5^0, equation ([59)1 is essentially the same to RPDE. 
Specifically, one only needs to set 

t/ = ^ exp {25U) , 

in order to transform equation ((59|) to ^{a, (3,U;n,m) = 0. For this reason, equation ([59]) will be 
referred to as RPDE for all values of parameter S. 

The function U may also be considered as a potential for ui and U2. To see this, we just have to 
solve S[U ; 6] for the derivatives of U. The compatibility of the resulting equations leads to the following 
system for ui, U2- 

2(ui — U2) / r2 / \ \ 

ui,a(3 = 7 To ToT ^ (wi Mi^ + 2d (n+ljmj 

2(U2 — Ml) / ^ r2 / 1 \\ 

M2,a/3 = 7 7^ TT7 ^ ("2,a M2,;3 + 2 d 71 (m + 1 ) 

+ (^(^,_,,)2_52(«_^)2 + ^) ((-+l)-2,„ + nu2,,) . (60ii) 
When S — Q, the last equations decouple easily yielding the RPDE pair 

£i{a, l3,ui;n + l,m) — , ^{a, (3, U2', n, m + 1) = . 

In case <5 ^ 0, system (pO)) may also be decoupled but the resulting equations are much more complicated. 
Specifically, one may solve equation (|60ip for U2 and substitute the result into equation (|60iip . This leads 
to a fourth order, second degree partial differential equation for ui, which is omitted here because of its 
length. Analogous considerations hold for the function U2- 

Remark 6.2 

System S[U;0] first appeared in T9j. A generalization of S'[[/;0] was derived in [531 [21] in the context 
of a symmetry reduction of the anti self dual Yang Mills equations. The relation of the latter to the 
Ernst- Weyl equation and the Painleve transcendents were also presented in [24j. □ 

Remark 6.3 

As already noted, S[U; 6] is integrable. A Lax pair for this system is given by 

^''^J—x[ U2 {m + U2U,) -U2 C/, J * " I -X)U, j * ' ^^^^^^ 



It can be obtained from the Lax pair (|26ip using the transformation ./# and performing the gauge 
transformation 

$ = (a - A)-"/2 _ A)-m/2 ^ 



□ 
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Remark 6.4 

It is worth mentioning that, equation ([59]) is the Euler - Lagrange equations 



da dp \dU^ap J da \dU^aJ 9(3 \dUjj 
corresponding to the Lagrangian 



For 5 = this reduces to the Lagrangian for the RPDE given in [12] . □ 
6.2 The system of H3 

The continuously invariant solutions of H3 are also related to solutions of RPDE. We establish this 
connection for the cases 5 — Q and 6 ^ Q separately. In each case, we introduce a potential function 
through a system of equations and use the latter to simplify the corresponding system Yi[u\. The 
resulting system can be decoupled leading to the RPDE for the potential function. 

6.2.1 Case I: (5 = 

Let us first introduce a potential ip for system corresponding to H3s=a- This is determined by the 
relations 

e^ in - aU^) e^ (m - pUp) 

V,a = ^ , ip.,p = z , (62) 

2ui 2u2 

where 

exp (— [/(a, /3)) — u{a,f3). 

We solve the above equations for ui, U2 and substitute the resulting expressions into the two first 
equations of Then, we perform the change of the independent variables 

(a, P) {a', P') 

and arrive at the following system for U and ■(/;: 

^'"^ - a ^ -p V^J' ^'^'^ 

2 

i',af3 = ^ _ p i.CtU^ai',/} - PU^f3lp,a) ■ (6311) 

On the other hand, using the above substitutions for u, ui and U2 and the change of the independent 
variables, the third equation of is identically satisfied by taking into account system (|63il) . 

The pair of equations (j63ip can be decoupled, and this leads to the following equation for the potential 

M{a, P,Tp;n,m) = 0. 

The decoupling can be achieved by solving equation (|63iip for one of the first order derivatives of U, 
e.g. U^cn and taking the compatibility condition between the resulting equation and (|63ip . The result 
is a relation for U^pjs. Finally, the compatibility condition between the latter and (|63ip implies that ip 
satisfies RPDE. 

On the other hand, system (|63ip may be decoupled leading to a fourth order, second degree partial 
differential equation for function U. First, we solve (|63ip for to get 
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where 



_ 2a/? (g - (3)U^0 + X 2{aU^ - (3AUp) 

a (m2 - /32 {/2j) ' a - P 



and 



The compatibihty condition dpip,a — 4'.ai3 imphes 

V',/3/3 = V',/? ( J 

FinaUy, the compatibility condition ^p.ap — da '>P,f3i3 leads to 

D„D^ In A = D„ (^^^ - 

If we write out the last equation explicitly, solve it for X and square the result, then we end up with 
a fourth order, second degree (in the highest derivative U.aapf}) partial differential equation. It is the 
modified partial differential equation (MPDE) presented in ,19j . 

6.2.2 Case II: 5 ^ 

In this case, we introduce the potential (f> by the relations 

^ Ui{nu + au,a) ^ U2{mu + (3u,p) 
a{uui A- 5 a) ' [3 {uu2 + 5 P) ' 

in view of which, ui, U2 may be expressed in terms of u, 4> and their derivatives: 

"1 \ 1 ' "2 = -5 ■ -, ^— ■ (64) 

We substitute the above relations into the two first equations of Yj[u] and set 
u{a,P) :— exp (— J7(a, /3) + (^^'(a,/?)) , 

(65) 

0(a, P) :~ —U{a,P) — S'ip{a, P) + n In a + m In /3 . 
Finally, the change of the independent variables 

(a, P) {a^ p^) (66) 

leads to the following system: 

U 1 ( 4a^U%-n^ t, 4/32t/2^ - m2 ^ \ ^ 

'''"'^ ^ i(^( — ^ — P — 

2 

V'.Q/S = {aU,atp,i3 - PU^pi;^a) ■ (67ii) 

a — p 

Moreover, the third equation of S[u] is satisfied identically, in view of transformations and by 

taking into account system (|67ip . 

We may decouple system (|67ip following the procedure described in the previous subsection. In this 
fashion, we conclude that, ip satisfies RPDE, as well. 
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Remark 6.5 

A Lax pair for system (|67i|) is given by 



1 / aU^a + A^JV.a 2 A -0,0 \ 



(68) 



The above equations follow from (|26ip by making the transformations (|64p - (|66p and, subsequently, 
performing the gauge transformation 



exp(^) 



□ 



6.3 Connection with previous results 

The preceding analysis shows that the integrable lattice equations H1-H3 and Ql are closely related, 
i.e. their continuously invariant solutions may be expressed in terms of solutions of RPDE. 

The relation between Hl^ H3s=o and Qls=o and RPDE was presented by Nijhoff, Hone and Joshi in 
[19j . starting from a different point of view. Specifically, the authors presented systems of differential- 
difference equations compatible with the above lattice equations, which, from our point of view, are the 
invariant surface conditions (|21iil21iii[) . Using the differential-difference equations, they constructed 
compatible systems of partial differential equations, which in turn lead to RPDE and MPDE. Actually, 
the systems appearing in the analysis of Nijhoff, Hone and Joshi do not differ from what we called 

To clarify this correspondence further, let us point out that, 

i) The system presented in |19| in relation with HI is actually S{u, — ui, —U2', 0). 

ii) In relation with Q1s=q, the authors of [19j presented a system of differential equations A, which 
can be decoupled leading to a fourth order partial differential equation, called Schwarzian partial 
differential equations (SPDE). Here, we have presented the corresponding system i.e. system 
(|48| . which may be decoupled for each involved function leading to RPDE. 

However, system A and the resulting SPDE are related to P5|) and RPDE, respectively. Indeed, 
starting from system (j48p . we make the change of the dependent variables 

2au a 2 j3u 3 

Ui ^ U + — — ^— - , U2 = U ' 



n (1 — ui) ' m (1 — U2) 

and, consequently, the change of the independent variables 



(a, P) 



1 1 

a ' l3 



This procedure leads to system A. Moreover, RPDE is mapped to SPDE using the above transfor- 
mation of the independent variables a, (3. 

iii) Finally, the authors of also presented the MPDE in relation with H3s=o, and a Miura trans- 
formation relating MPDE to RPDE. From our point of view, this Miura transformation is system 
(j63ip . which is equivalent to corresponding to HSs^q. 
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7 Conclusions and perspectives 



We have presented symmetry reductions of the Adler, Bobenko, Suris equations using both of the ex- 
tended three point generahzed symmetries admitted by them. Such reductions lead to special similarity 
solutions, which we named continuously invariant solutions. It was proven that these are determined 
by a system of partial differential equations, which is integrable in the sense that, it admits an 

auto-Backlund transformation and a Lax pair. 

The symmetry analysis and the corresponding reductions of system S [u] associated with the discrete 
potential and Schwarzian KdV equations led to new interesting results. In particular, it was shown that, 
the continuously invariant solutions of HI are determined by solutions of the continuous Painleve V and 
VI equations. Similar results and considerations were also presented with regard to equation Qlg^Q. 

We were also able to reveal the connection of T,[u\ to generating equations. In particular, we derived 
the generating equations to the H1~H3 and Ql members of the ABS family. In addition, we showed 
that, the continuously invariant solutions of HI, H3s=o and Ql s=o are related to RPDE, in accordance 
with the results of Nijhoff, Hone and Joshi in [19 . 

The construction of continuously invariant solutions of the other members of the ABS family and, 
especially, of the master equation Q4 is one of the interesting directions in which the present work can 
be extended. In addition, more general lattice systems possessing the consistency property can also be 
analyzed in the framework of continuously invariant solutions and generating equations. The discrete 
Boussinesq equation [2F and the discrete modified Boussinesq equation [26] are among the better known 
systems which can be brought into the above framework. 
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A The characteristic polynomials of the ABS equations 



HI : 


f{u,x,a) 


= 1 k{a,P)=(3-a 


G{x,y) = (x-yf 


H2 : 


f{u,x,a) 


^2{u + x + a) k{a,l3)=i3-a 


G{x,y) = {x-yf - {a 


H3 : 


f{u,x,a) 


= ux + ad k{a, P) = — (3 


2 G{x,y) = {ya- xP){yf3 


Ql : 


f{u,x,a) 


= {{u-xf -a^6'^)/a, k{a,(3) 


= -a(3{a - P) 




G{x,y) = 


a(3 {{x-yf -{a -(3)^6^) 




Q2: 


f{u,x,a) 


= {{u - x)'^ - 2a^{u + x)+ a^)loL 


, k{a, (3) = -a(3{a - P) 




G{x,y) = 


a/5((a;-y)2-2(a-/5)2(a; + y) + 


(a-Pf) 


Q3: 


f{u,x,a) 


4a(j2_i)(4a(au x){oix u) 


[a^ - 1)2^2) 




fc(a,/3) = 








G{x,y) = 


'"''4^''^ (4a/3(a^-/3x)(/3y- 


ax) + [a^ - p^)S^) 
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Q4 '■ f{u, X, a) = ((ua; + a{u + x) + .92/4)^ — {u + x + a){4:aux — g^)) /a 

_ ab{a^b+ab^ + [l2af3^-g2{a+2l3)-3g3]a+[l2f3a^~g2{0+2a)~3g3]b) 
K(a,P) - 4(q-/3) 

G{x,y) = {aoxy + ai{x + y) + a2){a2xy + a^ix + y) + ai) 

— {aixy + a2y + a2X + a'i)[aixy + d2X + a2y + a^) 

B Proof of Proposition 13.21 

Let u be a continuously invariant solution of the integrable lattice equation 

Q(w(o,o):"(i,o):'"(o,i)'"(i,i);"'/') = (69) 

and u be constructed in terms of u via the auto-Backlund transformation Bd(u,{t, A). The function u 
is another continuously invariant solution of equation (I69p , provided that it satisfies the system 

Q{u{o,Q)i'U'(^i,o),U{Q.i),U(i.i):,a,l3) = 0, (70i) 

— + ni?(u(o,o),W(i,o),U(-i,o),a) = 0, (70ii) 
d uiQ 0) 

— ^ "ii?({t(o,o),"(o,i)''"(o = 0. (70iii) 

Obviously, the first equation holds, since u is constructed using ]Bd(ii, u. A). It remains to show that 
equations (|70iill . (|70iiip also hold. 

To prove that (|70ii|l holds, we differentiate the first equation of Jid{u, u, A) with respect to a. Then, 
we use equation (|2Iiip and its shift in the n direction to substitute 9qU(q g) and daU(^i Q), respec- 
tively. Moreover, we use the determining equation for the generator Vi to substitute the derivative 
of Q(it(o,o) I W(i,o)! ^(o,o)i "^(1,0); ctj with respect to a. In terms of the above substitutions, we come up 
with the following equation 

{Q-M(o.o) + Q,iI(i,o) + "--R(u(0,0),U(l,0),'«(-l,0),a)^ = 0, 

where we have omitted the arguments of the function Q(u(o,o); W(i^o), U(o,o)j 'S(i.o); '^)- 
Finally, we eliminate the value U(i g) from the above equation and the latter becomes 

(G(u(o,0),M(l,0)) + (U(0,0),'«(0,0)) ^n) (^^(")~"^^ + ?^ ^ (-"(CO) , "(1,0) , "(-1,0) , ^ = 0. 

This equation involves the values of the function u and the value W(o.o) through the polynomials h, 
G. Thus, the corresponding coefficients of the various powers of U(o,o) must be identically zero. The 
matrix of the resulting algebraic system has rank two [13] , and the system admits only the zero solution, 
i.e. 

d u 



In the same fashion, we differentiate the second equation of the auto-Backlund transformation with 
respect to /3 and use equation (|21iii|) and the determining equation for the symmetry generator V2 to 
get that u also satisfies 

r{P) + "^ ^("(0,0), ■"(0,1), ■"(0,-1),/?) = 0. (71ii) 



Thus, function u satisfies system (|70ip . i.e. it is another continuously invariant solution. 
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C Proof of Proposition 13.31 

Proposition 13.21 implies that, if u is a continuously invariant solution, then function u, determined by 
]B(j(w, u, A), will be another solution of the same kind, and conversely. In other words, the functions u, 
u satisfy systems (|21l) and (j70ip . respectively. Thus, we can express the derivatives of ■U(o,o) in terms of 
the values and the corresponding derivatives of the function u and conversely. 
To achieve this, we solve 

Q(w(-i,o),W(o,o)7'"(-i,o),'«(o,o);a, A) = (72) 

for U(_i_o) E^nd (I70iip for U[-i.o)- In terms of the above substitutions, the fraction q) — U(_i.o)) 

becomes 

(r{a) , /,«(i.o)("(o.o),'"(i,o),a)\ Q,"(i,o) , Q:«(i.o)«(i,o) 
CaW(o,o) 7, 71 TP V ^ TT ' 

where we have omitted the arguments of Q(u(o.o)i '^^(1,0): ^'(o,o)j ""(i.o)! A). 

The derivatives of Q involved in the above expression are determined by the relations 

Q,u^i,o) _ fc(Q, A)/('U(o,o),M(i^o),Q;) 
Q,M(1,0) G(M(i^o): "(0,0): A) 

Q,u(i,o)U(i.o) _ 1 / k{a, A)/,ji(^_p) ("(0,0) . "(1,0) 7 a) ^ <^,"(i.o) ("(i.o): "(o.o) . A) \ 

Q,ii(i,o) 2 V G'(M(i,o),u(o,o),a, A) G(u(i^o), "(0,0), a, A) / 

which hold in view of the equation Q(u(o,o); ^^(1,0)5 ^*(o,o)7 'i*(i,o); Q!, A) — 0. 

We arrive at equation (|25iiip by substituting the above expressions into piiip and solving the resulting 
equation for 9q.U(o,o)- 

Conversely, we solve ((72)) for U(_i,o) and (j21iip for g). Then, we substitute the resulting ex- 
pressions into (|70iip and solve this equation for 9qU(o.o)- The final result is identical to the equation 
obtained by interchanging u and u in (j25iiip . 

Equation (|25ivp . i.e. the fourth equation of ©^(u, w. A), can be derived in a similar manner using 
equations (5(m(o,_i),M(o,o)> ""(0, -1)1 "(0,0); A) = and (I70iiip . 

Since the class of continuously invariant solutions is closed under and the initial solution u 
satisfies the same holds for the function u, i.e. the latter satisfies E[it]. Thus, Bc(w, u. A) defines 

an auto-Backlund transformation of system T.[u\. □ 
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